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Abstract 

The perturbative QCD approach to multiparticle production pre- 
dicts a characteristic suppression of particle multipUcity in a heavy 
quark jet as compared to a Ught quark jet. In the Modified Lead- 
ing Logarithmic Approximation (MLLA) the multiphcity difference 
dqe between heavy and hght quark jets is derived in terms of a few 
other experimentally measured quantities. The earlier prediction for 
5-quarks needs revision in the light of new experimental results and 
the improvement in the understanding of the experimental data. We 
now find 5be = 4.4 it 0.4. The updated MLLA results on 6be and 6c£ 
are compared with the present data from e~^e~ annihilation. Their 
expected energy independence is confirmed within the energy range 
between 29 and 200 GeV; the absolute values are now in better agree- 
ment with experiment than in the previous analysis, and the remaining 
difference can be attributed largely to next-to-MLLA contributions, 
an important subset of which are identified and evaluated. 



1 Introduction 



Since the early days of QCD, heavy quark physics has been one of the primary 
testing grounds for many aspects of the theory. In the last years a wealth of 
new important results on the profile of jets initiated by heavy quarks Q{b, c) 
has been reported by the experimental collaborations at LEP, SLC, Tevatron 
and HERA. Future progress is expected from the measurements at the LHC 
and a future linear e+e~ collider. These studies are important for the tests 
of the basic concepts of the QCD description of multiparticle production and 
also for the studies of new physics. 

Multiple hadron production in hard processes is derived from the QCD 
parton cascade processes which are dominated by gluon bremsstrahlung. An 
essential difference in the structure of the energetic heavy and light quark jets 
(£ = q = u,d, s) results from the dynamical restriction on the phase space of 
primary gluon radiation in the heavy quark case: the gluon radiation off an 
energetic quark Q with mass M and energy Eq ^ M is suppressed inside the 
forward angular cone with an opening angle ©o — M/Eq, the so-called dead 
cone phenomenon [1,2]. This is in close analogy with QED where the photon 
radiation is suppressed at small angles with respect to a primary charged 
massive particle. The suppression of the energetic gluon emission at low 
momentum transfer k± results, in turn, in the decrease of the heavy quark 
energy losses. This provides a pQCD explanation of the leading particle 
effect [3, 4] which is clearly seen experimentally in the bb and cc events in 
e'^e~ annihilation [5]; for recent reviews, see [6]. 

For a long time, there has been no clear explicit experimental visualisa- 
tion of the dead cone. Only recently, preliminary DELPHI results have been 
reported [7] which show the expected depletion of small angle particle pro- 
duction in 6-jets with respect to the heavy hadron direction, where the decay 
vertices of the heavy hadrons were reconstructed. Further detailed studies of 
the dead cone effect in different processes are needed. Some new results may 
come from the current analysis of the structure of the c-quark jets, produced 
in the photon gluon fusion in Deep Inelastic Scattering at HERA. 

It is worthwhile to mention that the difference in the radiation from mas- 
sive and massless quarks should also manifest itself in the QCD medium 
via suppression of the medium-induced radiative energy loss of heavy quarks 
propagating in a strongly interacting matter, see, for example, [8-11] and 
references therein. 

Studies of heavy quark jets are also important in the investigation of the 
properties of known or new heavy objects. For example, a detailed knowl- 
edge of the 6-jet profile is needed for the analysis of the final state in the tt 
production processes. Various aspects of studying new physics, in particular. 
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of the structure of the Higgs sector at the LHC and at a future hnear coUider, 
would benefit from the detailed understanding of the 6-initiated jets, see for 
example [12, 13]. 

The dead cone phenomenon leads to essential differences in the profiles 
of the light- and hcavy-quark-initiated jets. According to the concept of 
"Local Parton Hadron Duality" (LPHD) [14], the dead cone suppression of 
gluon radiation should result in the characteristic differences in 'companion' 
spectra and multiphcities of primary hght hadrons in these jets [1,2, 15]. 

In particular, as a direct consequence of the LPHD scenario, the difference 
of companion multiplicities of light hadrons in the heavy quark and light- 
quark jets at the same jet energy E'jet should be energy independent (up to 
a power correction O (M^/Ej^gJ), i.e. in e+e~ annihilation at c.m.s. energy 
W — 2£'jet one obtains the QCD prediction [6, 15] 

N^,{W) - N^qq{W) = const(W-). (1) 

The corresponding constant is different for c- and 6-quarks and depends on 
the type of light hadrons h under study. 

This prediction is in marked contrast with the expectation of the so-called 
naive model [16], which relates the multiplicities in light and heavy quark 
events based on the idea of the reduction of the energy scale, 

(^q) = ^^' ^-{xQ) = 0{as{W)). (2) 

In this case the difference of q- and Q-induced multiplicities would grow 
gradually with W as 

Ni,m - Ni^qiW) (X V^ln - • N^^^{W). (3) 

In this paper we focus on the analysis of the current experimental status 
of the difference of the average charged multiplicities 5u of events containing 
h- and light quarks in e^e~ annihilation in the available energy range. The 
situation with charmed quarks is considered as well. The main emphasis is 
on the comparison between the reanalysed data and the expectations based 
on the MLLA, in an extension of previous analyses [6, 15]. In addition, we 
discuss the size of the next-to-MLLA contributions. 

2 Theoretical Analysis 

Within the LPHD framework, the multiplicity of light hadrons in e"'"e~ anni- 
hilation events is proportional to that of bremsstrahlung partons. To predict 
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the QCD yield of light particles accompanying QQ production we have first 
to address the question on how the development of the parton cascade initi- 
ated by a heavy quark Q depends on the quark mass M. 



2.1 Structure of QCD Cascades in e^e QQ + . . . 

As well known, in the case of a light quark jet the structure of the parton 
branching of the primary gluon gi with energy uji (energy spectra, multiplic- 
ities of secondaries) is determined by the parameter 

Kq = Aujlsui^^, (4) 

where Gi is the angle between the gluon and the energetic quark. (The 
expression (jH) is written in such a way as to account for the next-to-leading 
correction due to large-angle soft gluon emission, up to the full jet opening 
angle 61 = vr, see for example [17,18].) For 9i -C 1 this parameter reduces 
to the gluon transverse momentum, ~ (c^iOi)^. The appearance of this 
scale is a consequence of colour coherence in multiplication of soft gluons 
which dominate the QCD cascades. This destructive coherence results in the 
Angular Ordering (AO) of successive parton branchings [19]. 

The corresponding parameter for a jet initiated by a heavy quark with 
energy Eq and the mass M reads 



2sin^) 



eo . ^. (5) 



Note that the same quantity hq determines the scale of the running coupling 
in the gluon emission off the massive quark. ^ 

The modification of the angular parameter in (jSJ caused by the heavy 
quark mass has a transparent physical interpretation.^ Consider radiation of 
a secondary gluon g2 with energy UJ2 <C oJi at angle Q21 relative to the primary 
gluon gi. Normally, in the "disordered" angular kinematics, 02i > Oi, the 
destructive interference between the emission amplitudes of g2 off the quark 
and gi cancels the independent radiation gi^g2 thus enforcing 

621 < 61 . (6) 



detailed analysis of the running coupling argument in the massive quark case can 
be found in Appendices to [3], see also [6,20]. 

^This argument is based on the discussion of two of the authors (YLD,VAK) with 
S.I. Troyan in the early 90's. 
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Meantime, in the massive quark case the interference contribution enters the 
game only when the angle B2 of (72 with respect to the quark is larger than 
the dead cone, O2 > ©o- Therefore, the cancellation leading to the AO 
condition (jH]) does not occur when the gluon gi is radiated inside the dead 
cone, 9i < 0o, and the jet evolution parameter © freezes in the 61 — 
limit. 

In physical terms what happens is the loss of coherence between Q and gi 
as emitters of the soft gluon (72 due to accumulated longitudinal separation 
Az > x\^^ ^ uj2^ between the massive and massless charges [vq ~ 1 — 6o/2 < 
1, Vi = 1). Indeed, during the formation time of the secondary radiation, 
t^P ~ {uj2&li)~^, the two sources — the quark and the gluon gi — separate 
in the longitudinal direction by 

Az ~ tf\vQ-ccose,\ ^ A[f)-%^. (7) 

It is the last factor that determines whether an interference is essential or 
not. When this ratio is larger than 1, the quark Q and gluon gi are separated 
enough for (72 to be able to resolve the two emitters as independent colour 
charges. In these circumstances gi acts as an independent source of the 
next generation bremsstrahlung quanta. Otherwise, no additional particles 
triggered by gi emerge on top of the yield determined by the quark charge 
(which equals the total colour charge of the Q + gi system) . 

In the massless quark case (Oq = 0) this consideration reproduces the 
standard AO prescription (jHl). In the massive quark case, the separation 
(incoherence) condition 621 < (Oi + Oq) results in (0) as the proper evolution 
parameter for the gluon subjet. 

The modification (0) may look superfluous since the soft gluon radiation 
inside the dead cone, 0i <C Oq, is suppressed. In spite of this, it is essential 
for keeping track of the next-to-leading order (MLLA) corrections in accom- 
panying multiplicities. In the Appendix A we recall the structure of the exact 
matrix element for gluon radiation off a heavy QQ pair and show how the 
parameter (jS)) naturally appears in the problem. 

2.2 MLLA prediction for accompanying multiplicity 
and its accuracy 

The light charged hadron multiplicity in heavy quark e~^e~ annihilation 
events at c.m.s. energy W can be represented as 

N^QiW) ^ iV,t_QQ (W) = iW) + <, (8) 
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where Nq is the charged muhiphcity of e^e events containing a heavy 
quark Q; [W) is the charged muhiphcity of hght hadrons accompany- 
ing the heavy quark production (excluding decay products of Q-flavoured 
hadrons) and Uq stands for the constant charged decay muhiphcity of the 
two leading heavy hadrons (n^"^ = 11.0 ± 0.2 for 6-quarks, n^'^ = 5.2 ± 0.3 
for c-quarks, see for example [15] for details of a previous evaluation). As 
shown in Appendix A, at W = 2Eq ^ M ^ Aqcd the companion multi- 
plicity Nqq{W) can be related to the particle yield in the light quark events 
e~^e~ qq {q = u, d, s) as [6, 15] 

N,,iW) - NQQiW) = N,,iV~eM) ■ [1 + O (a.(M)) ] , (9) 

where we approximately expressed the difference between the light- and 
heavy-quark generated multiplicities in terms of the light-quark event multi- 
plicity at reduced {W independent) c.m.s. energy Wq = \/eM, e = exp(l). 

Concerning the accuracy of Q, there are two separate issues one has to 
address, namely: 

1. the accuracy of the statement of the constancy of the l.h.s. of Q, 

2. to which accuracy this difference can be quantitatively predicted by 
means of pQCD (the r.h.s.). 

Left-Hand-Side. Answering the first question, it turns out to be insuffi- 
cient to compare particle multiplicities in a given order of perturbation the- 
ory. Indeed, within the next-to-leading accuracy (MLLA), for example, one 
takes into consideration ("exponentiated") y/a^ + ag effects in the anomalous 
dimension describing parton cascading, and 1 + a/oI terms in the normalisa- 
tion (coefficient functions). This allows to predict the l.h.s. of up to the 
NNLO correction the absolute magnitude of which is of the order of 

(Lh.s.) - (l.h.s.)*^^^^ = O {as{W) ■ N,g{W)) . (10) 

Meantime, the steep growth with energy of the multiplicity factor N{W) 
(faster than any power of In 14^) makes the neglected asN{W) correction 
dominate over the (presumably) finite r.h.s. in (jH)), thus endangering the 
very possibility of discriminating between Q- and g-jet multiplicities. 

However, examining the origin of perturbative corrections proportional 
to N{W) in (fTn|) one can see that all of them prove to be independent of 
the quark mass M, being inherent to the light quark jet evolution itself. 
For example, the first corrections of the order of as{W)N[W) to the MLLA 
expression (fTUI) come either from further improvement of the description of 
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the anomalous dimension A'~f{as) ~ Os^^ determining intrajet cascades, or 
from 0{asiW)) terms in the coefficient function due to 

• three-jet configuration quark + antiquark + hard gluon at large angle, 

• the so-called "dipole correction" to the AO scheme quark + antiquark 
+ two soft gluons at large emission angles (see [20]), 

both of which are insensitive to the 60 value with power accuracy O (Gq) <^ 1. 

In fact, the statement that the l.h.s. of Q does not depend on the an- 
nihilation energy follows from general considerations and should hold in all 
orders in perturbation theory with power accuracy, 1 + O {M'^/W'^). 

This is a very powerful statement which goes beyond the standard renor- 
malisation group (RG) wisdom about separation of two parametrically dif- 
ferent scales, W and M. Indeed, by looking upon the particle multiplicity 
as a moment (A^ = 0) of the inclusive fragmentation function, and by draw- 
ing an analogy with the OPE analysis of DIS structure functions (space-like 
parton distributions), one could expect for light- and heavy-quark initiated 
multiplicities 

^gl = f(M) = co..t(W) , (11) 

that is that their ratio rather than the difference is ly-independent. The 
RG motivated expectation (|TT| would have been correct if the quark mass 
M played the role of the initial condition for parton evolution — the trans- 
verse momentum cut-off. This is true enough for hard gluons with energies 
X = 2uj/W ~ 1 for which the region < M is indeed suppressed as com- 
pared to the massless quark case. It is not hard gluons that dominate the 
accompanying multiplicity however. 

Turning to (primary) gluons with x <^ 1 we observe that the radiation off 
light and heavy quarks remains the same down to much smaller transverse 
momentum scales namely, 

2M 

> uj- — = xM <^ M , 

which is nothing but the statement of the "dead cone" suppression discussed 
above. It is important to stress that, being based on the analysis of the ffist 
order gluon radiation matrix element, this conclusion is exact and holds in 
all orders in perturbative expansion. This follows from the fact that emission 
of gluons with x -C 1 is governed by the Low-Burnett-Kroll theorem [21] 
concerning the classical nature of soft accompanying radiation (following the 
dx ■ {Ijx — 1) distribution), which holds to power accuracy, see also [22]. 
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So the QCD coherence plays a fundamental role in establishing this re- 
sult [23]. Since the gluon bremsstrahlung off massive and massless quarks 
differs only at parametrically small angles G < Go, the AO (QCD coherence) 
then ensures that the accompanying cascading effects are limited from above 
by a finite factor N{W ■ Go) — N{M). A rigorous proof of the statement 
that W dependent corrections to the r.h.s of are power suppressed as 
M'^/W'^ (of subleading twist nature, in the OPE nomenclature) is lacking at 
the moment. 

By replacing the approximate MLLA multiplicities in pOj) by the exper- 
imentally observable multiplicities in it becomes possible to establish a 
phenomenological relation between the light and heavy quark jets with con- 
trollable accuracy. 

Thus, the difference in the mean charged multiplicities, 5q£, between 
heavy and light quark events at fixed annihilation energy W depends only 
on the heavy quark mass M and remains W^-independent (with power accu- 
racy) [15,23] 

<5q, = N^^ (W) - N^^ (W) = const (W^), (12) 
= A/'f (ly) - iVf (W^) = const (ly), (13) 

with Q = b,c and i = q = u,d, s. 

Right-Hand-Side. The r.h.s. of © is estimated with the MLLA accuracy. 
In general, this constant difference is proportional to N{M) and can be given 
in terms of the series in ^y as{M) as the pQCD expansion parameter. Let 
us remark that such an expansion formally relies upon treating the quark 
mass M as the second hard scale, as{M) <^ 1, and is bound to be only 
moderately satisfactory at best, since in practice, in (jH)), the bottom quark 
mass translates into ~ 8 GeV and the charm quark mass into Wq ~ 2.5 
GeV only.^ 

2.3 Quark mass effects in three-jet events 

Another powerful, and phenomenologically interesting, consequence of QCD 
coherence is that the structure of particle cascades in three-jet QQg events 
(with a hard gluon radiated at large angle) must be identical to that in 
the light-quark case everywhere, apart from the two narrow angular regions 
corresponding to the dead cones of the Q-quarks. More specifically, the 

■^Short-lived top quarks do not follow this pattern in the first place; the N^i notion 
being elusive, see for example [1,24,25]. 
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particle multiplicity in 3-jet events can be written in MLLA as the sum of 
quark and gluon jet multiplicities [26] (see also [27]) 

N,,,{W) = N,,{2E;) + \n,,{pI), (14) 

where E* denotes the g or g energy and p\ the gluon transverse momentum, 
both in the c.m.s. frame of the qq pair. Then, with Wqq = 2E*, we obtain 

^QQ.(^) - N^^-si^) = Nqq(^Qq) - NUW,,). (15) 

This may provide another handle for the detailed studies of the dead cone 
phenomenon at the reduced effective c.m.s. energies Wqq. 

2.4 Discussion and estimate of next-to-MLLA terms 
of the order of as{M)N{M) 

In [28] the equation ^ was evaluated exactly which constituted an attempt 
to improve the pQCD prediction beyond the accuracy beyond which 
(jni) does not actually hold. However, the main assumption of [28] that the 
companion multiplicity is generated by a single cascading gluon is not valid 
at this level. 

Next-to-MLLA correction terms are copious and it is hard to collect them 
all. There are, however, some specific contributions that look enhanced and 
can be listed and estimated. These are contributions containing an additional 
(semi- dimensional) factor vr^. 

In particular, to predict the event multiplicity at the agN level, one has 
to take into consideration large angle two soft gluon systems (aforementioned 
dipole configurations). This problem is discussed in Appendix A where we 
show how a vr^ enhanced correction emerges. Another correction of similar 
nature comes from the 1 — z rescaling of the argument of the dead cone sub- 
traction. This contribution is also extracted and analysed in Appendix A. 
It turns out to be numerically larger than the "dipole" contribution. These 
enhanced next-to-MLLA effects work in the same direction: they all tend 
to increase the difference between the light and heavy quark-initiated mul- 
tiplicities in 0. 
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3 Theoretical predictions confronted with ex- 
periment 



3.1 Experimental results on heavy quark multiplicities 
in e+e annihilation 

The experimental measurements of hadron multiplicities in bb and cc events 
produced in e~^e~ annihilation were performed in the wide range of c.m.s. 
energies ^/s = W from PEP, at = 29 GeV, to LEP2 at = 206 GeV 
[29-45]. For reviews on this topic see, for instance, [46-48]. Within the ex- 
perimental uncertainties the data clearly show that the differences 6i,e and 
6ce are fairly independent of the c.m.s. energy, as expected from the pertur- 
bative analysis, and in a marked contrast to the steeply rising total mean 
multiplicity N^'^^. 

This can be seen for example in Fig. 1, which shows a compilation of 
direct measurements of 6m, Eq. (fT^ . This figure is taken from [44] with the 
addition of the result from the VENUS experiment at ^/s = 58 GeV [36] as 
well as the preliminary result from DELPHI at ^/s = 206 GeV [49] . The 
dash-dotted line shown in Fig. 1 corresponds to the weighted average among 
all published results, Sl^^ = 3.12 ±0.14, assuming that the measurements are 
uncorrelated. 

It is worthwhile to mention that the first preliminary data on the multi- 
plicity difference between the 6-quark and udsc-quark large-angle 3-jet events 
produced in decays, are reported by DELPHI [50,51]. According to (fTHjl 
these results can be related to the multiplicity difference, 6bg', between the 
6-quark and the g'-quark {q' = u, d, s, c) events in e~^e~ annihilation mea- 
sured in the effective energy range Wqq ~ 53 — 59 GeV. The data points do 
not show any sizeable energy dependence and are consistent with the precise 
direct result from the VENUS experiment [36] at i/i = 58 GeV. 

As it can be seen in Fig. 1, within the experimental uncertainties most 
data points are consistent with the original MLLA prediction [15], S^'"^^ = 
5.5 ± 0.8. However, the precise results from the OPAL, SLD, DELPHI and 
VENUS experiments, which dominate the weighted average value 5^^^, are 
definitely lower. The "naive model", based on the reduction of the energy 
scale ^/s, which predicts the growing difference as in Eq. ^ and, therefore, 
the gradually decreasing 6be is strongly disfavoured. 



9 



DELCO/Mark II/TPC • OPAL 
□ TASSO * DELPHI 



* TOPAZ A SLD 

T VENUS O Markn 




L_i I I L_i I I I ■ ' i_J i_i I I I i_i L_i i_] I ■ ' I I i_i I I I i_i I I ' ■ I 

20 40 60 80 100 120 140 160 180 200 220 

Vs (GeV) 



Figure 1: Experimental measurements of 6m plotted as a function of the 
c.m.s. energy [29-45]. The 1992 MLLA expectation 6^^^^ = 5.5 ± 0.8 [15] 
(shaded area) includes experimental errors on nf"^ and light quark multiplic- 
ities at ^/s ~ 8 GeV. The prediction of the "naive model" [16] based on the 
reduction of the energy scale is also shown (dashed area). 



3.2 Test of MLLA predictions for b quark jets 

Our main goal here is to explain why the previous numerical value of the 
MLLA prediction of 6^^^^^ = 5.5 ± 0.8 [15] needs a revision. This value 
relies strongly on experimentally measured quantities, and some new relevant 
results became available since the analysis presented in [15]. Furthermore, 
we reanalysed the old data on charged multiplicities at low energies, which in 
addition to some small errors propagated in the literature until now affected 
the result presented in [15]. 

As we already mentioned, the difference between the MLLA result and the 
experimental data on 6be would allow to probe the size of the next-to-MLLA 
effects of order as{Mb)Ngg{Mb). 

Let us first take a fresh look at the MLLA expression for the charged 
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multiplicity difference 6m 



5; 



■MLLA 

'be 



nt-K^JV~eMb) 



(16) 



in order to establish whether and where the two terms in the r.h.s. of Eq. ()l(j|l 
require revision in the light of the current improvements in the understanding 
of experimental data. 

(i) Mean heavy hadron charged decay multiplicity, nf'^: 
In the analysis of Ref. [15] the average number of charged particles coming 
from the decay of two i?-hadrons was taken as nf^ = 2A/'f = 11.0 ± 0.2 
In the present analysis we used the most recent result obtained from the 
combination of the ALEPH, CDF, DELPHI, L3, OPAL and SLD data on B- 
hadron production [52], N^" = 4.955±0.062, with an addition of 0.485±0.065 
tracks to include the charged decay products of and A, as measured by 
OPAL [53]. There is still an issue of the role of heavier 5-hadron states 
{B*,B**,...) and on how fast the 'saturation' with rising energy can be 
established. Their contribution to the mean heavy hadron charged decay 
multiplicity, nf'^ is usually evaluated with the help of Monte Carlo models. 
We used the value 0.22 quoted by the SLD experiment [42], which should be 
almost independent on y/s for c.m.s. energies above the mass peak. 

We finally arrive at the value 



which practically coincides with the previous result in [15]. 

(ii) Subtraction term N^^{y/eMb): 
The second term in Eq. ()16|) is related to the radiation within the dead cone, 
where primary gluons emitted off the 6-quark and the 6-quark itself act as a 
source of secondary soft radiation. In order to quantify the size of this term 
we have to address first the issue of the definition of the 6-quark mass, which 
should be appropriate for the dead cone physics. 

As well known, within perturbative calculations it would be natural to 
take the pole in the quark propagator as the definition of the quark mass. By 
its very construction the pole mass is directly related to the concept of the 
free quark mass. However, due to the infrared effects the pole mass cannot be 
used with arbitrary high accuracy (see [54] for recent review and references). 
Though in a more sophisticated calculation a mass definition which is less 
sensitive to the small momenta may appear to be more appropriate, the 
uncertainties in the quark mass of order of Aqcd are far beyond the accuracy 
of our consideration here. So for the purposes of this paper we use the two- 
loop pole mass value, quoted in [5], 



dc 



11.10 ±0.18 



(17) 



{Mb)poie = 4.7 - 5.0 GeV 



(18) 
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which cover, in particular, some of the short distance mass prescriptions 
[55,56]. The scale = y/eMb at which the subtraction term N;jj{W^) must 
be evaluated is then ^ = (8.0 ± 0.25) GeV. 

Since there are no direct measurements of charged multiplicity at this 
energy, we estimate N^^ (8 GeV) in the following way: 

• use as many as possible experimental results on inclusive mean charged 
multiplicity iV^^d below and above ^/s = 8.0 GeV, rather than restrict- 
ing to a very limited energy range as in [15]; 

• fit the data points to evaluate Njl'^^{8 GeV) by interpolation, using dif- 
ferent parameterisations and over a wide energy range, in order to test 
the consistency and stability of the results and to estimate a reasonably 
conservative uncertainty for N^^di^ GeV); 

• evaluate N^^^ (8 GeV) by subtracting the c-quark contamination from 



We studied all available data on mean charged particle multiplicity, NfJ^^, 
collected in e'^e~ annihilations in the centre-of-mass energy range 1.4 GeV 
- 91 GeV. We considered only pubhshed results obtained in the continuum, 
thus away from the J/^ and T resonances, which were determined following 
what is now considered a standard convention [57], namely including in the 
evaluation of the mean value all charged particles produced in the decays 
of particles with lifetimes shorter than 3 • 10~^° sec. [58-74]. This means 
that the charged decay products of and of weakly decaying heavy-mesons 
(D,B, . . .) and baryons (A,E, . . .) as well as of their antiparticles must be 
considered, regardless of how far away from the interaction point the decay 
actually occurs. Unfortunately, some old publications, particularly those ob- 
tained at energies below 7 GeV, do not explain sufficiently well how the data 
were treated in this respect. We use only those ones which clearly considered 
at least charged decay products of Kg, that is known to be the dominant con- 
tribution at low energies. Furthermore, we do not consider results obtained 
at energies which might suffer from threshold effects due to charmed meson 
pair production, including higher mass states, notably the data collected by 
MARK I [59] in the interval 4.0 to 7.0 GeV. In order to evaluate N^adi^.O 
GeV) we fit the data points using the following parameterisations 



iV,t,(8 GeV). 




a + b-\n{s) + c-\n^{s), 



(19) 

(20) 
(21) 
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which are known to describe the data on mean charged multiphcity very 
well [75]. The parameters a, b and c, as well as the effective scale A not 
explicitly shown in (j2T| but that enters the definition of the running coupling 
as, are free parameters.'^ The two cases with three and five active flavours 
were considered in the calculation of as when making fits. 

We tested the consistency and the stability of the results by varying the fit 
energy range over the intervals: 7-14 GeV; 7-44 GeV (to include the results 
from PEP and PETRA); 7 - 62 GeV (to include results from TRISTAN) and 
7 - 91.2 GeV (to include results from LEPl and SLC), the common starting 
point of 7 GeV being well above the charmed meson production threshold. 
All mean multiplicities measured above ^/s = 10.5 GeV were corrected for 
the effects caused by the b-quark. At each energy the correction was ap- 
plied accounting for the fractions of the various quark species as predicted 
by the Standard Model and using the value 6m = 3.1 as measured experi- 
mentally. The total uncertainty associated with each data point was taken 
as the statistical and the systematic uncertainties added in quadrature. 

All fits give a very good x^, and the mean charged multiplicity predicted 
at i/i = 8 GeV is found to vary between 6.9 and 7.3. 

Our conclusion is that in the energy interval ^/s = 7.75 — 8.25 GeV 



The uncertainty includes the observed spread of values due to the choice of 
different parameterisations as well as the effect due to the uncertainty of the 
b-quark pole mass, {Mh)poie- 

The result presented in fl22|) . however, refers to a mixture of u,d,s,c- 
events, while for the determination of 5^1^^^ from Eqs. p6|] and (fTTj) only 
the contribution to from the light quarks (g = u,d,s), N^^, should 

be considered. In order to extract the light quark event multiplicity from 
(j221) we carefully studied the literature about the experimental results on the 
measurement of the multiplicity difference between the q- and c-quarks 



At the time of the analysis of Ref. [15] only the results from MARK II, 
TPC and TASSO were available. These results are affected by large un- 
certainties, and we also noticed in the literature some inconsistencies in the 
evaluation of S^^, which we corrected for. Much more precise results from 
OPAL [41] and SLD [42,45] are now available, and in the present analysis the 

^Note, that only (|21|l is pQCD motivated, but we are using (|19|l and (|2()|l as well for 
interpolation purposes and error evaluation. 



iV,td(8.0 GeV) = 7.1 ±0.3. 



(22) 



^ci c qq 



(23) 
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experimental value of 6^^^ to be used for the correction was reevaluated, as 
discussed in detail in Appendix B.l. It is shown there, that the experimental 
results from 29 GeV to 91 GeV are well consistent with a constant value, 
and a weighted average yields 



This value is about a factor two smaller than that used in [15], 6ce = 2.2 ±1.2, 
and is more precise. Since no direct measurements of at ^/s = 8 GeV exist, 
we assume its constancy also at lower energies, as in [15]. Clearly, a direct 
and accurate measurement of Sc£{^/s = 8) GeV, for example by analysing 
radiative events with the proper effective energy at the BaBar and Belle 
experiments, would be much desirable to validate our hypothesis. 

We finally correct iV^^^ for the effect of the 40% admixture of cc events 
using this new result on 5^^^, and find for the light quarks 



As a cross-check of this method, we estimate N^^ also in the following 
way. Besides the b-quark contribution, we subtract also the c-quark contri- 
bution from all the mean charged multiplicities measured above the c-quark 
threshold. This is done using the value of 6^^^ presented in ()24|1 and the 
Standard Model predictions for the c-quark fractions at each energy. 

With the exception of the MLLA parameterisation which in principle 
should not be used below the b-quark threshold, we can then extend the 
fitting procedure described above to the published results down to 1.4 GeV. 
This time the interpolation at 8 GeV provides directly the evaluation of N^^, 
to be used for the calculation of 5^1^^^. The values of Ng^{8 GeV) are found 
to range in the interval 6.45 - 6.65, completely consistent with the value of 
6.7 ±0.34 quoted in (j^. 

We also compared our findings with the results from several global QCD 
fits to N^ad 8.0 GeV. The numerical solution of the MLLA evolution 
equation for the particle multiplicity generated by light quarks, supplemented 
by the full 0{as) effects for e~^e~ annihilation [76], gives N^adi^-^ GeV) = 
6.5. In this fit no effort was undertaken to separate the contributions from 
different flavours, so the fit which includes low energy data as well should 
be placed in between Njl'^^ and A^^^, to be compared with (j22|l and (j25|) . 
The 3NL0-fit [77] using the data above 10 GeV gives N^adi^.O GeV) = 7.3, 
consistent with (j221)- Furthermore, a value A^^(8.0 GeV = 6.5 is found by 
running the Pythia 6.2 Monte Carlo program (in its default version) with 
light quarks only, with initial state radiation switched off and following the 




1.0 ±0.4. 



(24) 



N^liS.O GeV) = 6.7 ±0.34. 



(25) 
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standard convention for the definition of mean charged multiphcity,^ in good 
agreement with our result (j2Sl). 

Substituting (|T7jl and (j25|) into Eq. (|T6|l we arrive at the revised MLLA 
expectation for the multiphcity difference 

gMLLA ^ 4 4 ^ ^26) 

which is ~ 1.0 unit lower than the result reported in [15] and has half of its 
uncertainty. 

The comparison of the MLLA result ()26|) with the available experimental 
data on 6^ in e~^e~ annihilation is shown in Fig. 2; here we included also 
the reevaluated results of DELCO, MARK II, TPC, TASSO, TOPAZ and 
VENUS (see Appendix B.2). The new experimental average is given by 

= 3.14 ±0.14. (27) 

We can say that, qualitatively, the previous conclusion that the experimental 
mean value is lower than the absolute value of the MLLA prediction remains 
valid. Quantitatively, however, the agreement between the data and the 
theory definitely improves. 

Finally, we turn to the question of whether the remaining discrepancy 
can be attributed to the next-to-MLLA contributions. First we note that 
the experimental value of the multiplicity difference Nqq(W) — Nf^i(W) = 
nf^ — 51^^ = 7.96 ± 0.23 and the MLLA expectation Ngg{y/eM) as evaluated 
in ()25p differ by a relative amount < 20% which is of the order of the expected 
correction term in Q of (9 {as{M)). To gain insight at the quantitative level, 
we consider first the size of the above multiplicity difference in the Double 
Logarithmic Approximation (DLA). This is given by ()A.30|1 but with the 
r.h.s. replaced simply by Nqq[M) if the dominant contribution to Nq in Eq. 
()A.11|) is taken. For 6-quark jets, this requires evaluation of the multiplicity 
at Mb ~ 4.85 GeV. We estimate N^adi^b) ~ 5.1 (N^^^iMb) ^ 4.7) and, 
therefore, Sj^^^^ ~ 6.4. This is about two units above the MLLA prediction 
()26|1 which, in turn, is about one unit above the data in ()27|1 indicating 
convergence. 

In the next-to-MLLA two large "vr^-contributions" are derived explicitly, 
see Eq. f)A.30|) in Appendix A. The final expression involves the coupling 
at scale M which we derive from the 1-loop formula with A = 250 MeV, 
as typically used in MLLA applications (see for example, Ref [6,20]), and 
we obtain as{Mh) = 0.23 for Uf = 3 flavours. Then from ()A.30|1 we find 
Ng^iW) - N^QiW) = N^^{y/EM) X 1.27 ^ 8.5. This finally gives the result 

^T. Sjostrand, private communication. 
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Figure 2: Experimental measurements of 6be plotted as a function of the 
c.m.s. energy, ^/s■, data below 90 GeV reevaluated (see Appendix B.2). The 
revised MLLA expectation using 5^1^^^ = 4.4±0.4 is indicated by the shaded 
area. Also shown is the "naive model" [16] based on the reduction of energy 
scale (dashed area). 



including these next-to-MLLA contributions 6be ~ 2.6 ± 0.4. We, therefore, 
conclude, that the MLLA prediction is already close to the experimental 
data in ()27|) . and the remaining difference is of the order of the expected 
next-to-MLLA contributions. 

3.3 Results on charm quark jets 

Since the scale relevant for the charm quark, Wq ~ y/eMc, is significantly 
lower than in the 6-quark case the predictions are less reliable. 
The two-loop c-quark pole mass is quoted in [5] as 

(M,)p„,e = 1.47- 1.83 GeV. (28) 

We evaluated the size of the subtraction term N^^{Wq) where we followed 
the same strategy as described above, restricting the multiplicity fits to the 
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energy range 1.4 - 10.45 GeV. The predicted value at Wq = 2.7 GeV is found 
to vary between 3.5 and 3.9, and we arrive at N^^{2.7 GeV) = 3.7 ± 0.3. 
Using the c-quark decay multiplicity nf'^ = 5.2 ± 0.3 we obtain the MLLA 
expectation for the charged particle multiplicity difference in the c-quark case 

6^^^^ = 1.5 ± 0.4 (29) 

which is basically the same as the previous number S^^^"^ = 1.7±0.5 in [15]. 
The result fl^i is consistent with the new more precise experimental average 
given by Eq. (jSlJ. As in the case of 6be the theoretical MLLA result lies now 
above the experimental value which is expected due to the presence of the 
higher order effects. 

We also note an interesting aspect of the difference between the b and 
c quark multiplicities 6bc = Sm — Sd- Since the M-dependence of the next- 
to-MLLA term in ()A.30|1 is weaker than that of the leading Nqq{y/eM) con- 
tribution, the multiplicity difference 5cb = — See is less affected by this 
correction, and can be better approximated by the MLLA result. If we com- 
pare the experimental and theoretical numbers obtained from the results 
derived above 

^MLLA ^ 2.9 ± 0.6, C'' = 2.1 ± 0.4, (30) 

we find indeed that, contrary to the difference 6be, within the slightly larger 
errors, there is a reasonable agreement between the data and the MLLA 
prediction for this multiplicity difference involving b quarks. 



4 Conclusions 

The comparison of particle multiplicities in heavy and light-quark initiated 
jets provides a specific test of the perturbative approach to multiparticle 
production. In this approach the particle multiplicities in e~^e~ annihilation 
are directly proportional to the gluon multiplicities generated by multiple 
successive bremsstrahlung processes from the primary quarks. In the case of 
a primary heavy quark the small angle radiation is kinematically suppressed 
(dead cone effect). Also the subsequent gluon emission is affected by the 
mass effects, which results in the loss of coherence of soft gluon radiation off 
the heavy quark and the primary gluon. The result can be represented as 
an appropriate expansion in ^/a^ where the leading double logarithmic and 
next-to-leading (MLLA) terms have been known for quite a while, whereas 
certain large contributions in the next-to-MLLA order are discussed here. 

The main aim of this study is to sharpen the tests of the perturbative 
approach by accounting for all currently available data on e'^e~ annihilation. 
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More accurate theoretical predictions for the difference of multiphcities in 
hght and heavy quark jets are obtained. The expected energy independence 
of this difference is nicely confirmed. The same difference in 3-jet events is 
expected to agree with that in 2-jet events at the corresponding qq c.m.s. 
energy, and this is supported by preliminary data. As compared to the 
previous analysis, the updated MLLA prediction for the absolute value of 
the multiplicity difference comes closer to the experimental data. It is shown 
that the remaining difference is of the order of the next-to-MLLA corrections 
considered. This way the specific effects related to soft gluon bremsstrahlung 
from heavy quarks and their impact on the generation of the gluon cascade 
are reaching quantitative understanding within the perturbative approach. 

It would be very interesting to extend the measurements of 5u and 5ce. to 
lower energies, for example down to the region accessible at the B-factories. 
In particular, a direct measurement of Set. at a/s = 8.0 GeV, for example 
from the analysis of events with initial state photon radiation, would be very 
important to confirm our assumption that 5ce. remains constant below 29 
GeV. Further tests of the QCD predictions at higher energies at a future 
linear collider will be interesting. 
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Appendix A 



A.l Single gluon emission in MLLA and beyond 

The exact first order expression for probability of single gluon emission off 
the heavy quark pair can be written in the following form [3,22], in analogy 
with QED [78], 



Cpas dz dr] J ^ V - Vo , 2 
dwv = i 2(1 - z) — + z 



nv Z ^/l — T] I 7] 



^2 



1 1 

_r/ ^ 2 



Cv' , (A.l) 



with z the gluon energy fraction and r] an angular variable and 

Am? M 

l>7] = l-(3^ cos^ 0e > r/0 = YZ-^ ' ^ = ^ ^ ^ ' (^-2) 
where {3 is the quark velocity and 9c is the polar gluon angle in the QQ 



j3^ = l3\z) = l-^^ <v^ = l-Am^ > z. (A.3) 

The first term in curly brackets in ()A.1|) contains the main (double logarith- 
mic) contribution and corresponds to universal soft gluon bremsstrahlung. In 
accordance with the Low-Barnett-KroU theorem [21], both dz/ z and dz parts 
of the radiation density have a classical origin and are, therefore, universal, 
independent of the process (and of the quark spin). This term explicitly ex- 
hibits the dead cone phenomenon: "soft" radiation vanishes in the forward 
direction, sin 9c — *• 0, 77 — >• r]Q. 

The second term proportional to dz z (hard gluons) depends, generally 
speaking, on QQ production mechanism. Namely, both —1/2 subtraction 
term and the factor C,v = (3 — f ^) /2 = 1 + 2m? would be different for produc- 
tion current other than the vector current. We include this remark to stress 
that at the level of corrections (as well as of power suppressed effects 
O (yojm^)) the mean multiplicity acquires process-dependent contributions 
from 3-jet ensembles and cannot be treated any longer as an intrinsic char- 
acteristic of the QQ system. 

To obtain the mean parton multiplicity with the MLLA accuracy it suf- 
fices to supply ()A.1|) with the gluon cascading factor which depends, together 
with the running coupling, on the argument 

*?=(^V.,. (A.4) 
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see [3]. Neglecting relative corrections O {as) and O (m^) we obtain for the 
mean multiplicity 

Nqq{W) = No-Nu (A.5) 

where 



Cf <i2 A, 1 + (1 - ,.,,„, 

= ^ r'il /"*^2(l-2).|a,A/-dft). (A.7) 
We start by analysing the leading term ()A.6|) . 



Nq. The kinematical factor 1/^1 — rj somewhat enhances the contribution 
of the large angle region, rj = O {1), and should be taken into consideration 
in the leading DL term. The corresponding (SL) correction can be approxi- 
mately accounted for by pushing up the upper limit of the logarithmic inte- 
gration. Indeed, given that the factor JF = agN depends on 77 logarithmically, 
the chain of approximations follows: 

/' ^n,) = /' + /' ( 1 - 1) n,) 



t '^T{r]) + ln4 ■ ^(1) ^ —T{r]) . (A.8) 



It is straightforward to check that the omitted terms are limited from above 
by the O {as{W)) and O {m?) terms. Natural rescaling of the integration 
variable, t = rjW'^/i, leads to 

ATq = — / dz ^ / - [asUc] ih) , (A.9) 

with 

k'l = z^t. (A. 10) 

Now we represent f)A.9|) as 

iVo = N,^{W) - Ngg{M) - , (A.ll) 

where we have singled out an (enhanced) next-to-MLLA correction term that 
we will consider later, 

= N2{M) = ^ dz '^^' / ^ [aMa] (h) , (A. 12) 
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and introduced the function 

-1 ^ ^ ^^ ~^2 



= ^ / dz^-±^^^ [ ^ [a^ATc] {h) (A. 13) 

Jo ^ J t 

that describes the light quark event multiphcity at the c.m.s. energy W . 

We observe that the dead cone suppression naturally results in the expres- 
sion for the accompanying multiplicity in QQ events as a difference of light 
quark multiplicities at c.m.s. energies W and M. The MLLA correction A^^i 
defined in flA.7|) modifies the effective energy of the subtraction term N{M) 
in (jSH}. 



Ni. Since the r] integral in ()A.7|) is non- logarithmic and is concentrated in 
the region 77 ~ 770 <^ 1, we allowed ourselves to drop the l/-\/l — t] factor 
here as producing a negligible O {iv?) correction. We have 

A^i = — /' dz '^S^^ r ^ [aMc] (ho) , kl = -^M' u . (A.14) 
vr Jo z Ji l-z 

Though the collinear logarithmic enhancement disappears here, the soft one 
is still present (contrary to A^2) and promotes Ni to the a/oJ (MLLA) level. 

First we observe that the (1— -2) rescaling of the argument of the cascading 
factor [agAf] is negligible as it produces a O (^Q^s^^j correction (next-to-next- 
to-MLLA). Then, making use of the expansion 



r^Filnu) = F(0) + F'(0) + ... 



and replacing the factor 2(1 — z) by the numerator of the full quark gluon 
splitting function, 1 + (1 — 2;)^, we arrive at 

^1 = V {jy^ ^ '^\ asAfG]{zV^M) - ^[asAfcKM)^ , (A. 15) 

which holds with the next-to-MLLA accuracy (including 0{as))- By com- 
paring ()A.15|) with ()A.13|) and recalling ()A.10|) we can express the MLLA 
correction as the logarithmic derivative of the light quark multiplicity: 

Ni = 1 N'J^eM) ; N'^^iQ) ^ N,,iQ) . (A.16) 

Invoking ()A.5j) and (jA.llj) . for the QQ event multiplicity we finally obtain 



N,,{M) + ^N'^^{V^M) + 



2 



~ N,g{W) - N^^iV^M) - N2 . (A. 17) 

This proves the MLLA subtraction formula ©. 
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The effect due to the (1 — z) rescahng of the lower hmit of the t- 
integration in ()A.9|) produces a tt^ enhanced next-to-MLLA correction N2. 
We have 



-L / rf^ ^ L In [a^Uc] (M) 

TT Jo 2: l-Z 



TT V 3 4 



^""^ ^ ' [asMa] (M) • {1 + O {all\M))]. {K.li 



In terms of the event multiplicity {Nqq ~ ^CpjNc ■ No) we arrive at the 
relative correction 



N^iM) N^as /vr^ 5 



A. 2 Two gluon (dipole) correction 



(A.19) 



To derive the probabilistic MLLA equations describing parton cascades one 
has to analyse, in particular, ensembles of many energy ordered gluons radi- 
ated at arbitrary angles and demonstrate that, after having taken into full 
account multiple interference diagrams, one arrived at the pattern of angu- 
lar ordered (AO) successive gluon emission [20]. Reduction of interference 
graphs to the probabilistic AO scheme is not exact: there is a "remainder". 
In particular, the first such remainder appears at the a1 order and describes 
radiation of two soft gluons (with energies k2 <^ ki <^ W) at large angles 
with respect to the qq pair and to each other. The angular structure of the 
remainder R^"^"* is as follows 

= CpHl ■ N,Dl^^^^ + CpH^^ ■ N,Dl^_^^ , (A.20) 

where ± mark the momenta of q and g, the factor H describes independent 
gluon emission, 

HI = — , aik = g2 (Pt Pk) ^ ^ _ ^ 2 _ cos 0.^ ^ (A.21) 
aa [Pi ■ q) [Pk ■ Q) 

and D is the so-called "dipole factor" , 

^e[mn] = ^lm~^ln-: (A. 22) 

ji -\- Ctiffi fl£m / . r)n\ 

^im — ~~~~ • [A.Z6) 
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The dipole remainder possesses no collinear singularities, 

dilo f d^li 19 dx , ^, ^ vr^ , , ^ 

— - — -HiD^,,,.= nx = -C2= , A.24 

so that the integration of ()A.20|) over the gluon angles gives 

With account of the gluon cascading factor, logarithmic integrals over the 
gluon energies induce the next-to-MLLA correction to the event multiplicity 

Now we estimate the energy integrals using 



— ^^Gik') ~ -■^^G{k), 70 = w^^, 

k' 7o V TT 

with 7o the DLA multiplicity anomalous dimension, and obtain another tt^ 
enhanced relative correction 

NUQ) " 6 7o' 27r '24 ^ 

This means that the true multiplicity Ngg = Ng^^^^^ + ANgg and its MLLA 
estimate are related as follows 

^ N,M) . (1 + . g . (A.28) 

Now we return to the expression ()A.17jl . 

N,,iW) - NQQiW) = Nggi^eM) + N^{M) . (A.29) 

The first observations we make is that in the difference Ngq{W) — Nq-q{W) 
the two-gluon dipole corrections cancel since, as we discussed above, large 
angle soft gluon emission is insensitive to quark mass. Therefore, we can look 
upon the l.h.s. as being constructed of the true multiplicities. 

On the contrary, the factor Nqq on the r.h.s. of ()A.29|1 is the theoretical 
(MLLA) expression. Relating it with the true multiplicity via ()A.28|) results 
in 



Ngg{W) - Nqq{W) = Nggi^M) .1^1 + 



27r 



TT^ 5 



24 V 3 4 



(A.30) 
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where we inserted the expression flA.19|) for the first enhanced correction A^2- 
Numerically, the first term in the square bracket from the dipole corrections 
at large emission angles amounts only to about 4% whereas the second one, 
which improves the description of the small angle emission from the heavy 
quark, is about 5 times larger. The result ()A.30jl is not claimed to be complete 
at this order but it includes the important tt^ contributions considered to be 
dominant and shows the size of the next-to-MLLA terms. Remarkably, both 
corrections work in the same direction increasing the difference between the 
light and heavy quark companion multiplicities. 

Appendix B 

B.l On the measurement of 5ci 

The experimental determination of 5ci at different energies is very important 
for this analysis. As discussed in Sect. 3.2, a key point in the evaluation 
of the absolute value of the MLLA prediction for 5u is the determination 
of the light-quark mean multiplicity, N^^, at a/s = 8 GeV. Experimentally 
one measures the mean charged multiplicity of an unbiased inclusive sample 
of hadronic events, N^^, and then subtracts the contamination of heavy- 
quark-initiated events. This can be done if one knows the fractions of light 
and heavy-quark events present in the sample, fi and /q, as well as the 
difference between mean multiplicities of the heavy and light quarks, 5qi, 
using the relation 

Kaa = fi ■ + fc ■ (iV,t + ^cl) + h ■ (iV,t + ^^^)- 

At y/s = 8 GeV, where only the c-quark-initiated events are produced 
on top of the light-quark events, a direct measurement of 6ce is not available 
and, thus, its value must be evaluated from the knowledge of experimentally 
measured values of 6^ at different energies. Moreover, the knowledge of 6ce 
is necessary to derive Sm from the results of those experiments which do not 
measure directly the c-quark event mean multiplicity. The measurement of 
6ci is difficult because it is not easy to select experimentally a highly enriched 
sample of c-quark initiated events. 

So far, only five experiments published their results on the direct measure- 
ment of the mean charged particle multiplicities, N^^ and iV|'*, for e~^e~ —>■ cc 
and e+e^ il {ii = qq = uu, dd, ss) events, including the evaluation of sta- 
tistical and systematic uncertainties: MARKII [29] and TPC [32] at ^/s = 
29 GeV, TASSO [33,34] at = 35 GeV and OPAL [41] and SLD [42,45] 
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at ^/s = 91.2 GcV. These results, together with the derived values of Sd and 
their weighted average are presented in Table 1.^ 

It should be mentioned that the two results from SLD [42,45] were ob- 
tained from two completely independent event samples. The most recent 
one was collected with an upgraded detector, using a different experimental 
procedure and with different sources of systematic errors. We then consider 
the two results practically uncorrected . It should also be noticed that the 
results of MARKII and TPC presented in Table 1 are different from those 
derived in [15], and used to evaluate the light-quark charged mean multi- 
plicity at ^/s = 8 GeV in the same article. This is simply due to the fact 
that in [15] the values of N^^ used to calculate 6ce were not those quoted in 
the publications [29] and [32], but they were recalculated assuming a com- 
mon mean value for the total average multiplicity, N^^d^ as determined by 
different experiments at energies surrounding ^/s ~ 29 GeV. This procedure 
was meant to reduce the uncertainty on the derived values of Sd and 6be but 
is rather dangerous since information about the strong correlations among 
-^/lad' ^e^^ ^^"^ existing within the same measurement, is lost if 
one considers a mean value over different experiments for only one of these 
variables. That is why we rather preferred to use the published results which 
were all obtained within the same measurement. 

There are two more experimental results on N^^ and Nf^ published in 
the literature, one by the HRS collaboration at ^/s — 29 GeV [31] and 
one by the DELPHI collaboration at = 91 GeV [40]. Unfortunately, 
only the statistical uncertainties were evaluated in these analyses, and since 
the contribution of the systematic errors to the total error quoted by the 
other experiments is important, or even dominant, we did not consider the 
results from HRS and DELPHI in our weighted average. We show in the 
following that in any case, under reasonable assumptions about the size of 
the total errors, the final result would not change significantly if wc did. 
The DELPHI experiment measured N^^^ N^^ and Nf^ and found 5ci = 1.64. 
The total uncertainty (statistics and systematics combined) on 6be quoted in 
their analysis is about a factor two larger than those quoted by SLD [42,45] 
and OPAL [41], and if we assume a similar relative precision also for S^, by 
comparison with the SLD and OPAL total uncertainties we get for DELPHI 
See — 1-64 ± 1.2. Our weighted average in Table 1 would change to {6^) — 

^There are published results on the measurements of N^'^ and N^'^ also at LEP2 ener- 
gies [43,44,49]. Unfortunately, the limited statistics available at each energy did not allow 
the efficient c-quark tagging, comparable to that at the peak. Therefore, the selection 
of highly enriched c-quark samples was not possible. As a consequence the measurements 
of N^'^ are affected by large uncertainties, and cannot be used for a meaningful evaluation 
of 6ce- 
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Experiment 


x/i {GeV) 


pjch 
c 






MARKII [29] 


29 


13.2 ± 1.0 


12.2 ± 1.4 


1.0 ± 1.7 


TPC [32] 


29 


13.5 ± 0.9 


12.0 ± 0.9 


1.5 ± 1.3 


TASSO [33,34] 


35 


15.0 ± 1.2 


11.9 ± 1.2 


3.1 ± 1.6 


OPAL [41] 


91.2 


21.52 ± 0.62 


20.82 ± 0.44 


0.69 ± 0.62 


SLD [42] 


91.2 


21.28 d= 0.61 


20.21 ± 0.24 


1.07 ± 0.59 


SLD [45] 


91.2 


21.096 ± 0.653 


20.048 ± 0.316 


1.048 ± 0.718 


Average 








1.03 ± 0.34 



Table 1: Mean charged particle multiplicities, N^^ and Nf^, for cc and II 
{££ = uvl, dd, ss) events and the difference Sd = Nf'' — Nf^, measured at 
different energies. The results are corrected for detector effects as well as 
for initial state radiation effects. Charged decay products from Kg and A 
decays are included. We derived Nf^ for TASSO from the published values 
of N^^, N^^ and N^^^^ assuming the Standard Model quark fractions. The 
quoted errors are obtained by combining the statistical and the systematic 
errors in quadrature. OPAL and SLD errors on 5^ were published considering 
also correlations. The weighted average assumes no correlations among the 
various experimental results. 

1.07 ± 0.33 if we would include also this result. 

The HRS experiment measured Nf^ and Nf^, and found = 1.6. The 
results on Nf^ and Nf^ are consistent with those found by MARK 11 and 
TPC, and the size of the statistical errors are similar. If we attribute to the 
HRS value of 5^ a total uncertainty similar to those quoted by MARKII and 
TPC (here we assume a total error of ±1.5) and include also this measure- 
ment in our weighted average, we would get {S^i) = 1-09 ± 0.32. 

In conclusion, we use (5^) = 1.0±0.4 in the present analysis, and we point 
out that considering the current experimental precision there is no evidence 
of energy dependence of Sd in the range 29 GeV < ^/s < 91 GeV. 

B.2 About the measurement of 6bi 

In Table 2 we present an updated review of the experimental measurements 
of the mean charged particle multiplicities, Nj^^^, N^^ and Nf^, respectively 
for the inclusive sample (when measured), bh events and ll {ii — uu, dd, ss) 
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events. The difference 5^ = N^^ — Nf^ is also sliown. Tlie results are cor- 
rected for detector effects as well as for initial state radiation effects. Charged 
decay products from the Kg and A decays are included. The quoted errors 
are obtained by combining the statistical and the systematic uncertainties 
in quadrature. The published results on 5u from OPAL, SLD, DELPHI and 
VENUS take correlations into account. According to [46], the DELCO result 
appearing in table 2 was corrected by +25% as compared to the pubhshed 
DELCO data, i.e. 3.6 ± 1.5, to account for the overestimated h purity of the 
selected sample. 

We would like to stress at this point that the results on 5m presented in 
published compilations, including this one, are not all direct measurements. 
MARKII and TPC at = 29 GeV, TASSO at 35 GeV, OPAL, SLD and 
DELPHI at 91 GeV and DELPHI and OPAL at LEP2 energies, measured 
N^^^ N^^ and either Nfl^^, the inclusive mean charged multiplicity, or Nf^ (or 
both), from which 5u is calculated in a direct way. The other experiments, 
instead, have only measured N^^ and iV^^^, and, thus, one particular value 
for Nf^ or must be assumed in order to evaluate Nf- and 5u- In the 
previous reviews, the value of 5ct was the same as in [15], while in the recent 
publication by VENUS [36] the result of OPAL measurement [38] is taken. 
In Table 2 we used for all these experiments the new average value of 5ct 
presented in the previous section, 5ci = 1.0 ±0.4, and this explains why these 
results are not the same as those presented in previous publications. 
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Experiment 


V~s 


^had 






he 


DELCO [30] 


29 


12.3 ± 0.8 


15.2 ± 1.3 


11.6 ± 0.9 


4.5 ± 1.6 


MARKII [29] 


29 


12.9 ± 0.6 


16.1 ± 1.1 


12.2 ± 1.1 


3.9 ± 1.8 


TPC [32] 


29 




16.7 ± 1.0 


12.0 ± 0.9 


4.7 ± 1.4 


Average 


29 








4.4 ± 0.9 


TASSO [33] 


35 


13.4 ± 0.66 


16.0 ± 1.5 


11.9 ± 1.2 


4.1 ± 1.9 


TASSO [34] 


42.1 


14.9 ± 0.7 


17.0 ± 2.0 


14.3 ± 1.5 


2.7 ± 2.5 


TOPAZ [35] 


58 


14.21 ± 0.12 


16.24 ± 1.1 


13.57 ± 0.7 


2.7 ± 1.3 


VENITS [36] 


58 


16.79 ± 0.23 


19.38 ± 0.88 


16.07 ± 0.7 


3.31 ± 0.37 


MARKII [37] 


90.9 


20.9 ± 0.5 


23.1 ± 1.9 


20.1 ± 0.9 


3.0 ± 2.1 


DELPHI [40] 


91.2 




23.32 ± 0.51 


20.20 ± 0.45 


3.12 ± 0.68 


OPAL [41] 


91.2 




23.62 ± 0.48 


20.82 ± 0.44 


2.79 ± 0.30 


SLD [42] 


91.2 




23.14 ± 0.39 


20.21 ± 0.24 


2.93 ± 0.33 


SLD [45] 


91.2 




23.098 ± 0.378 


20.048 ± 0.316 


3.050 ± 0.311 


OPAL [44] 


130 




25.9 ± 1.3 


21.0 ± 1.4 


4.9 ± 1.5 


OPAL [44] 


136 




25.7 ± 1.7 


23.0 ± 1.6 


2.8 ± 2.0 


OPAL [44] 


161 




24.1 ± 1.7 


21.1 ± 2.1 


2.9 ± 2.3 


OPAL [44] 


172 




28.8 ± 2.2 


26.8 ± 2.1 


2.1 ± 2.5 


DELPHI [43] 


183 




29.79 ± 1.14 


25.25 ± 1.35 


4.55 ± 1.5 


OPAL [44] 


183 




28.3 ± 1.2 


26.8 ± 1.6 


1.5 ± 1.6 


DELPHI [43] 


189 




30.53 ± 0.78 


26.10 ± 0.97 


4.43 ± 1.05 


OPAL [44] 


189 




28.89 ± 0.77 


25.41 ± 1.0 


3.48 ± 1.2 


OPAL [44] 


192 




28.5 ± 1.4 


24.4 ± 1.9 


4.1 ± 2.0 


OPAL [44] 


196 




31.3 ± 1.5 


28.6 ± 1.6 


2.7 ± 1.7 


DELPHI [43] 


200 




29.38 ± 0.82 


25.99 ± 1.03 


3.39 ± 1.35 


OPAL [44] 


200 




30.3 ± 1.3 


25.6 ± 1.5 


4.7 ± 1.8 


OPAL [44] 


202 




29.9 ± 1.7 


25.5 ± 2.0 


4.4 ± 2.0 


DELPHI [49] 


206 




28.72 ± 0.77 


24.22 ± 1.09 


4.50 ± 1.17 


OPAL [44] 


206 




30.08 ± 1.0 


26.53 ± 1.4 


3.55 ± 1.2 
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Table 2: Corrected mean charged particle multiplicities and Sm at different 
energies (see text in Appendix B.2 for more details). According to [46], the 
DELCO result appearing in this table was corrected by +25% as compared 
to the pubhshed DELCO data, i.e. 3.6±1.5, to account for the overestimated 
b purity of the selected sample. 
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